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Abstract 

The multiplicity structure of the hadronic system X produced in deep-inelastic pro- 
cesses at HERA of the type ep — > eXY , where Y is a hadronic system with mass 
My < 1.6 GeV and where the squared momentum transfer at the pY vertex, t, 
is limited to \t\ < 1 GeV 2 , is studied as a function of the invariant mass Mx of 
the system X. Results are presented on multiplicity distributions and multiplic- 
ity moments, rapidity spectra and forward-backward correlations in the centre-of- 
mass system of X. The data are compared to results in e + e~ annihilation, fixed- 
target lepton-nucleon collisions, hadro-produced diffractive final states and to non- 
diffractive hadron-hadron collisions. The comparison suggests a production mech- 
anism of virtual photon dissociation which involves a mixture of partonic states 
and a significant gluon content. The data are well described by a model, based on 
a QCD-Regge analysis of the diffractive structure function, which assumes a large 
hard gluonic component of the colourless exchange at low Q 2 . A model with soft 
colour interactions is also successful. 
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1 Introduction 



The observation of "Large Rapidity Gap" (LRG) events JT],§] in deep-inelastic ep scattering 
(DIS) at HERA, which are mainly attributed to diffractive photon dissociation UK], has led to 
renewed interest in diffractive phenomena (5j and how they can be understood within quantum 
chromodynamics (QCD). This paper continues a series of final-state studies by HI and 
ZEUS |K)pfL"2]] and presents the first results on the charged particle multiplicity structure in 
LRG events. 

Following conventions adopted in earlier HI analyses [f7, 13 1, the large rapidity gap events 
studied here are experimentally defined in terms of the generic process ep — ► eXY, where the 
hadronic systems X and Y are separated by the largest rapidity gap in the event; Y is the system 
closest to the outgoing proton beam @,|,13|. Events with no activity in a large pseudo-rapidityQ 
domain adjacent to the outgoing proton beam are selected. In these events the proton either 
stays intact (Y = proton), or is excited to a low- mass system. The system Y has longitudinal 
momentum close to that of the proton beam and small transverse momentum. In such events, 
the system X, measured in the central part of the detector, can be viewed as mainly resulting 
from the dissociation of a photon with a virtuality Q 2 . The requirement of a large rapidity gap 
implies that the invariant masses Mx and My of the systems X and Y are small compared to 
W, the centre-of-mass energy of the "y*p system. In the following, the terms "LRG events" and 
"diffractive events" will be used synonymously, although, in practice, some contribution from 
non-diffractive and double-diffractive processes is to be expected, particularly at large Mx- 

In diffractive DIS it is useful to view the interaction in a Lorentz frame where the target 
proton is at rest. For small values of Bjorken-x, XBj, the virtual photon fluctuates far upstream of 
the proton target into a quark-antiquark pair which can subsequently evolve into a more complex 
partonic system (qq,qqg,- ■ •) before the actual interaction occurs p4-16|. The virtual photon 
can thus be described as a superposition of Fock-states with different partonic content [17]. 

Diffraction in general [18], and dissociation of the photon in particular, arises from the fact 
that the strength of absorption of the various Fock-states depends on the internal degrees of 
freedom and quark-gluon composition of the dissociating object [19|-21|. Specifically, whereas 
the total 7*j> cross section measures the average absorption strength, the magnitude of the 
diffractive cross section, on the other hand, is related to its fluctuations 1£,20|. The diffractively 
produced hadronic final states are therefore expected to carry information, not only on the 
parton composition of the virtual photon wavefunction components and their interaction with 
the target, but also on their respective contributions to the diffractive cross section. 

Virtual photon dissociation can also be approached from a t-channel perspective. In Regge 
phenomenology, the interaction takes place via a factorisable exchange of the pomeron (iP) and 
of reggeons related to mesons. It has been suggested to endow the pomeron (and reggeons) 
with a partonic sub-structure and to use the concept of parton distributions in the pomeron 
and in sub-leading reggeons to model diffractive deep-inelastic scattering P2]| . In the proton's 
infinite-momentum frame, the pomeron (or reggeon) has a fraction xjp = (q- (P — P'))/ (q ■ P) of 
the proton's four-momentum P (q and P' being the four-momentum of the exchanged photon 
and the system Y, respectively), while the fractional momentum of the exchange carried by the 
struck parton is /3 = XBj/xjp. This approach, adopted in [13] assuming the DGLAP evolution 
equations [23], indicates, as was already conjectured in 24], that the pomeron must have a 
large hard gluon content at low Q 2 but that a sub-leading meson exchange is also needed at 
larger values of xjp or Mx, if the basic hypothesis of factorisation of each component is to be 
maintained. 

The parton structure of the colourless exchange, as deduced from an analysis of the total 
diffractive DIS cross section, can be tested in studies of diffractive final states |8|,^,|l2, 25|. These 
confirm the need for a pomeron dominated by hard gluons at a starting scale Qq ~ 3 GeV 2 . 

1 H1 uses a laboratory coordinate system with the z-axis aligned with the proton beam direction. Pseudora- 
pidity is defined as r\ — ~ ln(tan |), where 9 is the polar angle with respect to the proton direction. 
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This paper complements previous work and presents an analysis of charged particle mul- 
tiplicity distributions, multiplicity moments, charged particle density in rapidity space and 
forward-backward multiplicity correlations measured in the centre-of-mass (CMS) of the sys- 
tem X. The emphasis is placed on a comparison with data from e + e~ annihilation, fixed-target 
DIS, hadron diffraction and soft non-diffractive hadron-hadron collisions. Monte Carlo models 
which represent various theoretical views on diffraction in DIS [26,27] are also confronted with 
the data. 



2 Experimental procedure 
2.1 The experiment 

The experiment was carried out with the HI detector |28| at the HERA storage ring at DESY. 
The data were collected during the 1994 running period when 27.5 GeV positrons collided with 
820 GeV protons, at a centre-of-mass energy of 300 GeV. The following briefly describes the 
detector components most relevant to this analysis. 

The energy of the scattered positron is measured with a backward electromagnetic lead- 
scintillator calorimeter (BEMC), extending over the polar angular range 151° < 6 < 177° with 
full azimuthal coverage. The BEMC electromagnetic energy resolution is &e/E ~ 0.10/ y/E\GeV] 
© 0.42/i?[GeV] © 0.03, while the BEMC energy scale for positrons is known to an accuracy of 
1% | |29| |. A backward proportional chamber (BPC), situated immediately in front of the BEMC 
and with an angular acceptance of 155.5° < < 174.5°, serves to measure the impact point of 
the scattered positron and to confirm that the particle entering the BEMC is charged. Using 
information from the BPC, the BEMC and the reconstructed event vertex, the polar angle of the 
scattered positron can be determined to better than 1 mrad. A scintillator hodoscope behind 
the BEMC is used to reject beam-induced background based on a time-of-flight measurement. 

The hadronic final state is measured by tracking detectors surrounded by calorimeters. The 
Central Tracker consists of inner and outer cylindrical jet chambers, z-drift chambers and 
proportional chambers. The jet chambers, mounted concentrically around the beam line, in- 
side a homogeneous magnetic field of 1.15 Tesla, provide both particle charge and momen- 
tum measurement from track curvature and cover an acceptance region defined by the an- 
gular interval 15° < 6 < 165° and transverse momentum px > 0.1 GeV. Up to 56 space 
points can be measured for tracks with sufficiently large px- The resolutions achieved are 
ct Pt /pt ~ 0.009 • px [GeV] © 0.015 and ag = 20 mrad [28, 3C] with a track finding efficiency 



above 95% for tracks well contained in both jet chambers. In addition, forward going particles 
can be detected by the Forward Tracker in the polar angular range 8° < < 20°. The liquid 
argon calorimeter (LAr) extends over the polar angular range 4° < < 154° with full azimuthal 
coverage. The LAr hadronic energy resolution is (Je/E ~ 0.50/ \J i?[GeV] © 0.02 as determined 



in test beams [31]. A study of the transverse momentum balance between the hadronic final 
state and the scattered positron has shown that the absolute hadronic energy scale is known 
to an accuracy of 4%. Backward going hadrons can be detected by the BEMC. The hadronic 



energy scale of the BEMC is known to a precision of 20% |29|. 

Forward energy at small angles is observed in several detectors near the outgoing proton 
beam direction. Particles reach these detectors both directly from the interaction point and as 
a result of secondary scattering with the beam pipe and other adjacent passive material. The 
detectors are thus sensitive to energetic particles produced in directions that are beyond their 
geometrical acceptances. The effective ranges of sensitivity to energy flow are 3.5 < i] < 5.5 for the 
copper /silicon sandwich (PLUG) calorimeter, 4.5 < ij < 6.5 for the Forward Muon Spectrometer 
and 6.0 < ij < 7.5 for the Proton Remnant Tagger, which consists of scintillation counters and 
is located 24 m from the interaction point |J. These detectors overlap considerably in their 
rapidity coverage, thereby allowing for intercalibration of their efficiencies. 
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2.2 Monte Carlo simulation of LRG events 



Monte Carlo generated LRG events obtained with the RAPGAP 2.02 generator |27j are used 
to correct the observed distributions for detector acceptance, inefficiencies and smearing effects. 
This model correctly accounts for many final-state features of LRG events (§|,|9, 12, 13], but was 
not tuned to the multiplicity data presented in this paper. All generated events go through a full 
simulation of the HI apparatus and are passed through the same analysis chain as the real data. 
The detector simulation is based on the GEANT program The RAPGAP predictions shown 
in subsequent figures were derived from a model version which incorporates the recent results 
of the QCD-Regge analysis of the diffractive structure function [13]. Results are also shown 
from the DIS Monte Carlo event generator LEPTO 6.5 |33| and from the JETSET parton- 
shower model [34] for simulation of e + e~ annihilation hadronic final states. Important technical 
aspects of these models are given in recent HI publications (8|,|9],[l3| and are not repeated here. 

The RAPGAP generator models LRG events as deep-inelastic scattering of a virtual photon 
off a pomeron or reggeon coupled to the initial-state proton. The pomeron and reggeon are 
given a partonic content. The meson structure function is taken to be that of the pion [35|. 

Different partonic sub-processes are implemented using Born-term and first-order pertur- 
bative QCD matrix elements: eq — ► eq scattering, QCD-Compton scattering (eq — ► eqg) and 
boson-gluon fusion (BGF) off a gluon in the colourless exchange (eg — * eqq). Their relative 
contributions are controlled by quark and gluon densities of the exchange as determined in 
the HI QCD-Regge analysis [13[ of the LRG event cross section using the DGLAP evolution 
equations. In this model a "pomeron remnant" heads in the direction opposite to the virtual 
photon, consisting of a quark or a gluon for eq or eg scattering, respectively. The fragmentation 
of the partonic systems created in the sub-processes eq — ► eq, eq — > eqg, is thus expected to be 
analogous to that in e + e~ annihilation with centre-of-mass energy ^/s = Mx- However, in BGF 
the initial partonic system consists of a pomeron remnant (gluon) and a qq pair in a colour-octet 
state. This process has no equivalent in e + e~ annihilation. 

To assess the sensitivity to the quark-gluon content of the pomeron, results are presented for 
two sets of parton distributions (labelled hereafter "RG F 2 D (fit 3)" and "RG (fit 1)"): i) a 
"hard gluon" distribution ( "fit 3" in |L|] ) whereby gluons carry > 80% of the momentum at the 
starting scale Qq = 3 GeV 2 ; ii) a "gq-only" distribution whereby only quarks are present at Qq 
("fit 1" in (l^]). The latter model version is disfavoured in fits to |R| and by diffractive 

final-state studies reported in |8l|s|,|l2|p|. 

Higher-order effects in the QCD cascade are treated with the parton shower model, as im- 
plemented in LEPTO [36|. Hadronisation is carried out with the Lund string fragmentation 
scheme, as in JETSET 7.4 [34]. QED radiative processes are included via an interface to the 
program HERACLES f37|. 

Soft colour interactions form the basis of an alternative model, which is implemented in the 



LEPTO 6.5 [26 1 generator. In this scheme |pq1 , LRG events are the result of a normal deep- 
inelastic scattering on the proton followed by a long-time-scale random colour rearrangement in 
the soft field of the target. This mechanism leads to a fraction of events with a large rapidity 
gap of the type studied in this paper. 

Studies of diffractively produced meson systems in meson-hadron interactions with masses 
Mx below 10 GeV have revealed striking similarities with e + e~ final states at yfs = Mx (see 
e.g. |39|| ). To analyse virtual photon dissociation along similar lines, the JETSET parton shower 
model for this process is used (labelled hereafter "JETSET-e + e~"). The predictions are cal- 
culated for a standard mixture of primary (u, d, s) quark-antiquark pairs only. A primary cc 
component (e + e~ — > cc) has been neglected, motivated by the prediction that, for a quark- 
dominated pomeron, heavy-quark production is suppressed |4(|. It was verified, however, that 
inclusion of a contribution from primary cc pairs has little effect on the multiplicity and would 



2 The predictions from fit 2 in [[13| are very similar to those of fit 3 for the observables studied here and not 
shown. 
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not alter the conclusions of this analysis. The generated e + e~ events are rotated such that the 
z-axis coincides with the initial qq axis of the event in the CMS. Rapidity^ is calculated rela- 
tive to that direction. The use of JETSET-e + e~ permits a consistent comparison with the HI 
data and avoids difficulties due to the differing experimental treatment of e.g. strange particle 
production in different e + e~ experiments. 

To study the systematic uncertainties arising from background processes, several other Monte 
Carlo generators are used. 



The DIFFVM generator [41] models the low mass region of diffraction (Mx < 1.1 GeV) 
via production of the vector mesons p(770), uj(782) and </>(1020). This model further includes a 
simulation of proton diffractive dissociation. 

The PHOJET generator is used to estimate background from photoproduction processes. 
The model simulates non-diffractive reactions, elastic vector meson production, vector meson 
production with proton dissociation, single-photon diffractive dissociation and double diffraction 



and is in broad agreement with experimental results at HERA [43|. 



2.3 Event and track selection 

A neutral current DIS event selection is made [44| by demanding a well-reconstructed scattered 
positron detected in the BEMC with an energy, E' e , larger than 10 GeV. A subsample of domi- 
nantly diffractive events is then selected by requiring no activity above noise levels in any of the 
forward detectors or in the most forward part (rj > 3.2) of the LAr Calorimeter. 

Further cuts are applied to ensure that a positron is detected and reconstructed with high 
quality. An event vertex, reconstructed from tracks in the central trackers, within z = ±30 
cm of the mean vertex position, is required to reject beam-induced background. Events with a 
time-of-flight veto from the scintillator hodoscope are rejected. 

The standard deep-inelastic kinematic variables (ysj, Q 2 ) are reconstructed with the meth- 
ods described in |7|,|l3| and which use both the scattered positron and the hadronic final state. 
The kinematic variables (xjp, Mx) characterising the final state in LRG events, are obtained 
from a combination of tracker and calorimeter information with an algorithm for track-cluster 
association which avoids double counting. In |7],[l3| it is demonstrated that Mx is adequately 
reconstructed across the kinematic range of the measurement with a resolution of about 25%. 

Corrected data are restricted to a kinematic region where the acceptance of the HI detector 
is high and the contribution from non-DIS background is low. Together with the requirement 
E' e > 10 GeV, an upper limit on ysj at 0.6 ensures that the photoproduction background is 
less than 0.3% in the selected event sample. The lower limit ysj > 0.05 ensures substantial 
hadronic energy flow in the detector and adequate resolution in jjBj- The ysj cuts correspond 
roughly to a range 70 < W < 230 GeV. Non-diffractive contributions are suppressed by requiring 
xjp < 0.05. A lower cut on Mx of 3 GeV excludes light vector meson contributions. In addition, 
the requirement of an absence of activity in the forward detectors imposes the approximate 
restrictions My < 1.6 GeV and \t\ ~ 1 GeV 2 , although these variables were not measured 
directly. The kinematic regions to which the data are corrected are summarised in Table [lj; 
there \t m i n \ is the minimal kinematically allowed value of \t\. The event sample consists of 4738 
events, corresponding to an integrated luminosity of 1.3 pb _1 . 

The multiplicity analysis is based on charged particles. Only tracks observed within the 
acceptance of the central tracking detector, and which are successfully fitted to the primary 
event vertex, contribute to the uncorrected multiplicity. Further details on the track selection 



criteria and efficiencies can be found in |44]. 



The multiplicity distribution is corrected with an iterative matrix migration method based 
on full Monte Carlo simulation of the detector response. The method is described in detail in [fi4| , 
45 1 . The results are cross-checked with a fit of a Negative Binomial distribution, smeared for 



3 Rapidity of a particle is defined as y — 0.5 In [{E + p z )/(E — p z )], where E is the energy and p z the momentum 
component along the direction of a predefined axis; the pion mass is assigned to each particle. 
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Quantity 


Lower limit 


Upper limit 


Q 


7.5 GeV 2 


100 GeV 2 


VBj 


0.05 


0.6 


X]p 


U.UUUo 


U.UO 


M x 


3 GeV 


36 GeV 


\t\ 


I train \ 


1 GeV 2 


My 


proton mass 


1.6 GeV 



Table 1: Limits of the kinematic regions considered for corrected data. 



detector acceptance and inefficiencies, to the observed multiplicity distribution. This parametric 
technique is known to be less sensitive to the generator input, but yields only the lowest moments 
of the multiplicity distribution. For the measurement of rapidity spectra, a standard bin-by-bin 
correction procedure is used, with cross-checks provided by the matrix-unfolding methods. 

Charged decay products of Kg, A and A and from weakly decaying particles with lifetimes 
larger than 8 • 10 -9 s are subtracted from the multiplicity distribution through the unfolding 



method [44]. Hadrons associated with the target remnant system Y are excluded from the 



multiplicity measurement. 
2.4 Systematic errors 

Several sources of possible systematic errors are investigated. The analysis is repeated for each 
source and the changes to the results are added in quadrature. For illustration, the typical 
systematic error on the mean total charged multiplicity and on the central rapidity (—0.5 < y < 
0.5) particle density are given in square brackets for each source separately. 

• The error due to the uncertainty of the energy scale of the hadronic final state is esti- 
mated by scaling the LAr, BEMC and Central Tracker energies by ±4%, ±20% and ±3% 
respectively [0.4%, 0.4%]. 

The systematic uncertainty in the reconstruction of the scattered positron is studied by 
varying the energy E' e and polar angle 6' by ±1% and ±1 mrad, respectively [0.3%, 0.9%]. 

• The influence of the pomeron (reggeon) flux and pomeron (reggeon) structure function 
used in the Monte Carlo generator for correction is investigated by reweighing the (3, xp 



and t distributions for Monte Carlo events as in [13] [f3: 1.3%, 2.2%; xjp: 0.3%, 1.6%; t: 
1.1%, 3.5%]. 



Both the colour dipole model, as in ARIADNE [46], and the parton shower model are used 
to evaluate the influence of these event generation schemes on the corrections. The full 
difference is taken as the systematic error [0.6%, 0.( 



The strangeness-suppression parameter (PARJ(2) in JETSET j34|), affecting the rate of 
strange particle production in the simulation of the fragmentation process, has been varied 
in the range 0.2-0.3 according to recent results on strange particle production |47],^8| [0.1%, 
0.1%]. 

Track-quality criteria (such as the track length and the number of hits) are varied to esti- 
mate systematic errors related to an imperfect description of the acceptance and efficiency 
of the Central Tracker in the Monte Carlo simulation [1.3%, 3.1%]. 

An uncertainty of 30% is assumed on the Monte Carlo correction outside the tracker 
acceptance (9 lab outside the range 8°-165° or p l f < 0.1 GeV; the range between 8°-15° 
has been cross-checked with data from the Forward Tracker) [3.6%, 11 
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• Background events are suppressed by the event selection criteria. Remaining background 
contamination is estimated by including events simulated with the PHO JET and DIFFVM 
generators in the Monte Carlo event sample [PHOJET: < 0.1%, < 0.1%; DIFFVM: 0.7%, 
1.0%]. 

The number of events with initial and final-state QED radiation is changed by ±50% [0.4%, 
0.1%]. 

• A fit to a smeared Negative Binomial distribution is used as a cross-check on the matrix- 
unfolding results [error on central rapidity particle density 2.8%]. 

3 Results 

All data presented belowQ are corrected for the effects of acceptance and resolution of the HI 
detector in the kinematic ranges specified in Table |l|. To optimize the statistical precision of 
the different measurements presented here, a fine- and coarse-grained binning in Mx is used. 
The requirement of a forward rapidity gap ensures that the hadronic final states of the system 
X are well contained in the central detectors. The data span the Mx range from 3 to 36 GeV, 
distributed over the intervals as listed in Table |2[ Statistical and systematic errors on the data 
points shown in the figures are combined in quadrature. Unless otherwise stated, where two 
error bars are displayed the inner one is the statistical and the outer shows the total error. In 
comparing HI LRG data at a given Mx to data from other processes, the corresponding centre- 
of-mass energy scale is chosen to be W for fixed-target DIS data, y/s for e + e~ and non-diffractive 
hadron collisions, Mx for hadro-produced diffractive states. 



Mx range 


(M x ) 


(P) 


(Q 2 ) 


no. of events 


(GeV) 


(GeV) 




GeV 2 




3-8 


5.4 


0.41 


21 


1492 


8-15 


11.4 


0.17 


26 


1515 


15-30 


21.1 


0.06 


27 


1359 


4-6 


5.0 


0.43 


22 


638 


6-8 


7.0 


0.30 


23 


530 


8-11 


9.5 


0.21 


26 


737 


11-15 


13.0 


0.13 


26 


778 


15-19 


16.9 


0.08 


27 


543 


19-24 


21.3 


0.06 


27 


468 


24-36 


29.1 


0.03 


27 


562 



Table 2: Corrected average Mx, (3 and Q 2 and the number of observed events for the different 
intervals in Mx considered in this analysis. 



3.1 Multiplicity moments 

The lowest-order moments of the multiplicity distribution, the average multiplicity (n), the 
dispersion D = ((n — (n)) 2 ) 1 ^ 2 and the normalised second-order factorial moment R2, have 
been measured as a function of Mx- The moment R2 is defined as R2 = R2/(n) 2 with 

R 2 = (n(n-l))=J2 p nn(n-l). (1) 

n 

4 All data are available in numerical form on request and can be retrieved from the Durham HEPDATA 
database. 
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The latter quantity is equal to the integral of the inclusive two-particle density over a given 
domain in phase space and is a measure of the strength of correlations among the produced 
hadrons p^ , |49[| . Note that in the case of uncorrelated particle production, the probability to 
produce n particles, P n , follows a Poisson distribution with the result that R2 = 1. 

Fig. |l|a shows the dependence of the mean charged particle multiplicity on Mx in full phase 
space. The HI LRG data can be parameterized by a form (n) = a\ + C12 In M x + 03 In 2 M x , with 
a x = 2.2 ±0.4, a 2 = 0.08 ±0.17 and a 3 = 0.21 ±0.02 ( X 2 per degree of freedom = 0.4; statistical 
errors only), indicating that (n) increases faster than the logarithm of the centre-of-mass energy 
of the system. In non-diffractive DIS at HERA, a similar rate of increase is observed with respect 
to W [ 44 1 . Also shown are measurements of (n) for the diffr actively produced system X in the 
reactions ir^p — > X^p and K + p — > X + p pS|, |5C[|. Although the two data sets agree well for 
Mx ~ 10 GeV, (n) in LRG events exceeds that in meson diffraction at larger masses. There are 
no meson diffraction results with Mx ~ 15 GeV. 

The meson-diffraction data are close to the e + e~ annihilation results, represented here by the 
predictions from the JETSET parton shower model (dotted line) which is known to reproduce 
well the e + e~ multiplicity data over a wide energy range (see e.g. [^l|,^2|). 

Results on the dispersion and the correlation parameter R2, (Figs. ||b,c) confirm that also the 
second-order moments in LRG data are similar to meson diffraction and e + e~ for Mx ~ 10 GeV 
within the precision of the measurements^. Stronger multiplicity fluctuations and correlations 
than in e + e~ are observed at larger Mx- The rise of R2 with Mx shows that KNO-scaling [ |53[ ] 
does not hold in the Mx range studied here. 

The similarities seen in Fig. [l]a-c have led to the view |3^, [54|] that in meson diffraction 
the (mainly longitudinal) momentum exchange with the target leads to an excited meson state 
which can be pictured as a colour-string of invariant mass Mx stretched between the valence 
q and q of the meson. This string subsequently hadronises in a way similar to a quark pair in 
e + e~ — > qq at the corresponding centre-of-mass energy = Mx- A comparative study of the 
thrust distribution, energy and quantum number flow in the rest-frame of the system X and 
also in e + e~ further support this interpretation ]pq| . For Mx ~ 10 GeV the same idea has been 
succesfully applied to proton diffractive dissociation assuming that the baryonic system X now 
results from the fragmentation of a (valence) quark-diquark string, thus explaining observed 



similarities with DIS lepton-nucleon data at values of W comparable to Mx [39, 59. Due 



to the larger values of Bjorken-x involved, the latter reaction is dominated by quark-diquark 
fragmentation. 

Combining these experimental results, it follows that low-mass diffraction (above the reso- 
nance region) in hadron collisions and photon dissociation in DIS at HERA may be interpreted 
as the hadronisation of a single string, or colour dipole, with colour triplet-antitriplet endpoints. 
The larger multiplicity moments seen in LRG events for Mx ~ 10 GeV relative to the other 
processes suggest, however, that the above interpretation is incomplete and that high-mass 
diffraction involves additional mechanisms. Indications from experiment on the possible nature 
of these mechanisms exist for hadron diffraction. 

In high-mass proton diffraction, measurements show that the multiplicity structure of the 
system X deviates from the expectations for quark-diquark fragmentation and becomes, in fact, 
similar to that of soft non-diffractive interactions with y/s = Mx [57|. 

Within the framework of the Dual Parton Model (DPM), which is phenomenologically very 
successful ] p8|] , soft non-diffractive collisions are described by the fragmentation of two or more 
strings corresponding to single or multiple pomeron exchange in the elastic channel. The sim- 
ilarity between non-diffractive and diffractive processes is explained in DPM |5£| by assuming 
that the colourless exchange in the latter becomes resolved in a qq pair at large Mx and sub- 
sequently interacts with the dissociating hadron. The diffractive state is then described by two 

5 The errors on the NA22 data points for R2 are derived from published results for (n) and (n(n — 1)) — (ri) 2 
whereby the (unknown) correlation between these quantities is neglected. The quoted errors are therefore to be 
considered as a lower limit on the true error. 
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colour strings, one stretched between a valence quark of the excited hadron and a quark in the 
exchange, the other between a diquark and the remaining quark of the qq pair. 

Multi-string systems are known to lead to a faster than logarithmic increase of (n) with 
energy, to wider multiplicity distributions and stronger long-range particle correlations than 
single-string fragmentation p8fl . It is also important to note that only colour triplet-antitriplet 
strings are considered in DPM. 

Additional mechanisms, besides qq fragmentation are included in present models for diffr ac- 
tive DIS. Fig. |l|a-c show model calculations with RAPGAP (fit 3) (solid line) which describe 
the data well. The difference between RAPGAP and JETSET e + e~ follows from the presence, 
in RAPGAP, of additional diagrams involving gluons from the colourless exchange, leading to a 
large contribition from boson-gluon fusion. The partonic state in lowest-order BGF consists of 
a gluon (the "pomeron remnant") and a qq pair in a colour-octet state. The fragmentation of 
this state allows for various string topologies, including two-string configurations, thus leading 
one to expect further similarities with large Mj hadron dissociation. 

The admixture of the BGF sub-process with the qq and QCD-Compton processes naturally 
explains larger mean multiplicity and stronger fluctuations. On the other hand, the results for 
RAPGAP (fit 1), with a quark-dominated pomeron leading dominantly to qq parton states, are 
very similar to those of JETSET e + e~, as expected, and are not shown. 

The data can also be qualitatively understood in the photon dissociation picture of diffrac- 
tion. The lowest-order ("aligned jet" ]i~5| |) excitation (7* — > qq) is dominant for M x < Q 2 
and leads to a final state similar to that in e + e~ annihilation. In addition, higher-order fluc- 
tuations (such as 7* — > qqg where the gluon has low momentum), which resemble the BGF 
sub-process, are believed to contribute at larger Mx and to effectively interact as an octet-octet 
colour dipole (l£ , 60] . Due to the octet colour charge at the dipole end-points, such a system 
hadronises with a larger mean multiplicity than a qq state [61]. 

The LEPTO model with "soft colour interactions" (dashed curve) is seen to also agree with 
the HI data although it predicts somewhat larger multiplicity fluctuations above Mx ~ 20 GeV. 
This model too contains a sizeable BGF contribution. However, diffraction is viewed here as a 
final-state interaction and does not invoke the notion of colour-neutral exchange. 

The moments of the multiplicity distribution for particles with positive and negative rapidity 
("forward" and "backward", respectively) are displayed in Figs. [j]d-f. Rapidity is calculated in 
the rest-frame of the system X (the "7*iP" centre-of-mass system) with the positive longitudinal 
momentum axis pointing in the 7* direction^, assigning the pion mass to each charged track. 

The HI data show no evidence for an asymmetry between forward and backward hemispheres, 
in contrast to what is observed for the mean multiplicity measured in fixed-target fip DIS 
(Q 2 > 4 GeV 2 ) [p2| , where the influence of proton fragmentation on the backward-hemisphere 
multiplicity distribution is known to be substantial. 

The fip data in the current fragmentation region, where the comparison with LRG data is 
most relevant, agree well with the predictions for e + e~, as expected for production dominated by 
quark jets. The LRG results are also here characterized by larger (n) and stronger fluctuations 
above M x Z 10 GeV. 

The RAPGAP and LEPTO models also predict forward-backward symmetry of the single- 
hemisphere moments and describe the data adequately. The earlier noted differences with e + e~ 
annihilation for the full phase space moments are also seen here. 

The role of gluons in high-mass photon dissociation is prominent in all the DIS models con- 
sidered, in contrast to models for hadronic soft diffraction where quark (diquark) fragmentation 
is dominant. The indication from the models that the pomeron-remnant has a large gluon con- 
tent also opens interesting opportunities for comparison with gluon-jet fragmentation in other 
processes. 

6 The pomeron direction cannot be unambiguously determined since the outgoing system Y is not measured. 
Since its transverse momentum is small, it has been assumed that the IP direction is collinear with the incident 
proton in the rest frame of the system X. 
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3.2 Multiplicity distributions 



The multiplicity distributions in full phase space have been measured, separately for negatively 
and positively charged tracks, in three intervals of Mj. The results are displayed in Fig. @a-c in 



the form of a KNO-distribution [53|: ip(z) = (n)P n plotted as a function of the normalised mul- 
tiplicity z = n/{n). No significant difference is observed between the distributions for positively 
and negatively charged hadrons. The data are well reproduced by the RAPGAP model (solid 
curves) although there are indications that it underestimates the high-multiplicity tail of the dis- 
tribution at large Mx- The comparison with JETSET e + e~ predictions (at yfs = (Mj)) shows 
that the multiplicity distribution is broader in the LRG data, indicative of stronger correlations 
among the hadrons. The predictions for LEPTO are similar to those of RAPGAP. 

Figs. |2|d-f further illustrate the forward-backward symmetry of the system X, now for the 
all-charged multiplicity distribution. The RAPGAP and LEPTO predictions are also forward- 
backward symmetric but tend to fall below the data at large z. The single-hemisphere dis- 
tributions are closer to the e + e~ expectations than in full phase space (cfr. Figs. §a-c). This 
difference can be understood as the effect of correlations between hadrons emitted in opposite 
hemispheres which, as will be shown below, are larger in the DIS LRG data. 

3.3 Rapidity spectra 

The charged particle rapidity density in three intervals of Mx is shown in Fig. ||[ The spectrum 
rises slowly with Mx in the central region and a rapidity plateau develops with increasing phase 



space. These features confirm earlier observations in hadron diffraction [55, 57, 63[ that the 
diffractive system hadronises in a jet-like manner both in the forward and backward regions M, 
H, O]. There is no evidence for a significant forward-backward asymmetry of the y-spectraQ 



contrary to what is observed in /j,N interactions [64, 3Jj|. 

The particle density in the central region is much larger in LRG events than in /j,N interac- 
tions for W values close to {Mx)- It is also larger than in e + e~ annihilation (at ^/s = (Mx)) 
according to the JETSET expectation. The RAPGAP (fit 1) model curve is close to JETSET 
and the fixed-target data and predicts a particle density which is too low. Both RAPGAP with 
a hard gluon distribution and LEPTO describe the rapidity spectra, although there are small 
deviations in the lowest Mx bin. 

Fig. ^a further compares the Mx dependence of the central particle density (defined as the 
mean multiplicity in the region —0.5 < y < 0.5) in LRG events to e + e~ expectations, to that 
in iiN collisions |pi,|65|, non-diffractive meson-proton |3(| collisions and proton diffraction [67] . 
The particle density near y = is seen to be larger in LRG events than in all the other processes. 

The excess particle production relative to that in e + e~ and fiN indicates that additional 
mechanisms besides hard and soft gluon bremsstrahlung from quarks are needed (cfr. Sect. 3.1). 

The comparison with non-diffractive meson-proton and high-mass proton diffraction further 
shows that the central particle density in processes which are believed [58 to involve two or 



more strings with colour triplet-antitriplet end-points (qq and quark-diquark strings) is also 
significantly lower than in the LRG data. This, together with previous observations, argues 
in favour of models which attribute a higher gluonic content to the partonic system created in 
virtual photon dissociation than in the other processes. 

An estimate of the importance of an additional gluonic component may be obtained by 
assuming that the particle density in the central region is a linear superposition of two contribu- 
tions, one arising from qq fragmentation (including additional QCD radiation), the second from 
a colour octet-octet string configuration. This hypothesis is in line with expectations from the 
photon dissociation picture of diffractive DIS (see e.g. p3[). 



7 The rapidity spectra have also been recomputed with rapidity defined along the thrust axis in the CMS of 
the system X (not shown). Except for the lowest Mx interval, no significant difference is seen with the results in 
Fig. g| This is consistent with the observation in || that the thrust axis in LRG events is strongly aligned with 
the 7* direction in the 7*iP system. 
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Using the EMC data for the former, and JETSET simulations of a colour-singlet gluon-gluon 
string for the latter, it is found that, at (Mx) = 11-4 GeV, approximately equal contributions 
of the two components are needed to explain the particle density at mid-rapidity|^|. 

The RAPGAP (fit 3) predictions for the central particle density are shown in Fig. |I|a. They 
are compatible with the LRG data only above Mx > 10 GeV and are nearly Mx independent. 
The LEPTO model, on the other hand, predicts a rather stronger dependence on Mx, closer to 
the tendency observed in the HI data. The enhanced particle density, both in RAPGAP (fit 3) 
and in LEPTO, are related to the large contribution from boson- gluon fusion. The RAPGAP 
(fit 1) model version follows closely the JETSET e + e~ prediction. 

In order to investigate the sensitivity of the results to possible contributions from non- 
diffractive processes, the analysis has been repeated changing the cut xp < 0.05 to xp < 0.025. 
Within errors, no significant effect on the results was observed. 



3.4 Forward-Backward Correlations 

In this section, differences between LRG events and final states in other processes are further 
examined through a measurement of the correlation between hadrons emitted in opposite event 
hemispheres. These so-called "forward-backward" correlations are known to be sensitive to finer 
details of the fragmentation process and, in particular, to the presence in an inclusive event 
sample of several distinct sub-classes of events |68j . 

In previous experiments, the forward-backward correlation was analysed by studying the 
regression between the forward multiplicity, rip, and the backward multiplicity, ng- The corre- 
lation is usually well parameterised by a simple linear dependence 



(rip) = a + b- n B . (2) 

Forward-backward correlations have not previously been measured in diffractively produced 
final states. For reasons of statistics, matrix techniques as used in other analyses |[l], [52| ^] 
to unfold the two-dimensional forward and backward multiplicity distributions have not been 
employed. Instead, the forward-backward correlation parameter is estimated from the separately 
unfolded and corrected multiplicity distributions in full phase space, in the forward and in the 
backward hemispheres. Exploiting the relation between the dispersion for the full phase space 
(D) and that for the forward and backward hemispheres (Dp and Db), one can define the 
correlation parameter p as 

D 2 - Dl -Dl 
-= 2D F D B ■ < 3 > 

The parameter p is identical to the slope b in eqn. (|2[) in the case of forward-backward symmetric 
systems |7Cfl . 

Fig. |3|b shows the parameter p in three intervals of Mx for the LRG data. Also shown are 



data on the parameter b for p,p collisions [62], for non-diffr active tr^/K^p collisions compiled 
in [70 1 and JETSET predictions for e + e~. 

In spite of the large errors, there is clear evidence for stronger correlation in LRG events 
than observed in e + e~ annihilation and in the p,p data for energies above > 10 GeV. At 
lower energy, phase space effects are important and mask possible differences in dynamics. The 
correlation strength in diffractive DIS is comparable to that in meson-proton interactions. 

At LEP, where a value of 6 ~ 0.1 is measured, OPAL finds that the small correlation 
observed in an inclusive sample of e + e~ events is primarily due to the superposition of events 
with distinct numbers of jets and, therefore, different average charged multiplicity. Sub-classes 



of n-jet events (n > 2) show no or even negative correlations [51, 52|. In vp and vp reactions 



8 This is consistent with the contribution of about 50% from boson-gluon fusion events to the total diffractive 
cross section estimated with the RAPGAP model; the latter value depends, however, on the pt 2 cut-off value, 
here chosen to be 2 GeV 2 [0 . 
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no clear evidence for correlations is observed. These data therefore show that forward-backward 
correlations are small at energies where the production mechanism is believed to be dominated 
by single-string qq or quark-diquark fragmentation. 

In contrast, abundant data from hadron interactions, compiled in [70|, which cover the range 



10 < yfs < 900 GeV, show that the correlation increases logarithmically with energy, with b 
as large as 0.65 ± 0.01 at yfs = 900 GeV. The strength and energy dependence of the effect is 
attributed to strong event-to-event fluctuations of the particle density as occur e.g. in the multi- 
string Dual-Parton model due to fluctuations in the number of strings and strings overlapping 



in phase space |p8 | 



The observation of forward-backward correlations in LRG events with a strength comparable 
to that in soft hadron interactions adds further support to the view that the inclusive sample 
of DIS LRG events is a mixture of states with distinct hadronisation properties. To disentangle 
their precise nature and relative contribution, more differential studies will be needed, however. 

In present models for diffractive DIS, distinct production processes are readily identified 
and related, either to the differences in parton composition and absorption probability of vir- 
tual photon Fock-states, or to hard quark- and gluon- initiated interactions off a colourless ex- 
change. That a mixture of such contributions leads to significant forward-backward correlations 
is demonstrated by the predictions for the parameter p from RAPGAP-i 7 ^ (fit 3) (solid line) 
and LEPTO (dashed) which are close to the HI data for Mx ~ 10 GeV. The large difference 
between these DIS models and JETSET for e + e~ illustrates the sensitivity of this correlation 
measure to differences in the dynamics of these two processes. 

4 Summary and conclusions 

The charged-particle multiplicity structure of large-rapidity-gap events of the type j*p — > XY 
in deep-inelastic scattering at HERA has been measured. The major fraction of these events 
is generally interpreted as due to diffractive dissociation of the virtual photon on the proton, 
7*p — > Xp. 

Multiplicity distributions, lower-order moments, rapidity spectra and correlations between 
hadrons emitted in opposite hemispheres in the rest-frame of the system X have been presented 
as a function of the invariant mass Mx ■ 

The data have been compared with e + e~ annihilation (at yfs = Mx), lepton-nucleon data 
inalf range comparable to the M^-range in the HI data, with hadro-produced diffractive final 
states, and also with data from non-diffractive hadron-hadron collisions at yfs ~ Mx- The main 
observations are the following. 

• The mean total charged particle multiplicity (n) is a function of Mx and increases propor- 
tionally to In 2 Mx- The inclusive rapidity spectrum is forward-backward symmetric in the 
rest-frame of A. A plateau develops with increasing Mx- Both (n) (for Mx ~ 10 GeV) 
and the particle density near y = are larger than in DIS at comparable values of W, than 
in e + e~ annihilation at yfs = Mx and than in hadro-produced diffractive final states. Fur- 
thermore, the particle density in this central region is also higher than in non-diffractive 
collisions at yfs = Mx- 

• For Mx ~ 10 GeV, multiplicity fluctuations are larger than in e + e~ annihilation and than 
in the current fragmentation region of lepton-nucleon interactions at comparable values of 
yfs and W, respectively. The forward-backward multiplicity correlations are also larger 
and of comparable strength to those measured in hadron interactions at yfs = Mx- 

The distinctive characteristics of large-rapidity-gap events mentioned can be globally under- 
stood if it is assumed that the photon dissociation mechanism involves a mixture of different 
partonic states wherein gluons play an increasingly important role as Mx increases. A large con- 
tribution from gluon-rich states is also required to explain the steep rise with increasing 1/xsj, 
at fixed Q 2 , of the diffractive as well as the total virtual-photon proton cross section [16,72|. 
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Good agreement with the data is achieved with a model which assumes that the diffractive 
process is initiated by the interaction of a point-like virtual photon with a gluon-dominated 
colour-singlet object emitted from the proton, as is suggested by a perturbative QCD-Regge 
analysis based on DGLAP evolution of the diffractive structure function. However, some devi- 
ations are seen in the large- n tail of the multiplicity distribution at high Mx ■ 

A model with soft colour interactions which rearrange the colour topology after a normal 
deep-inelastic scattering also describes the data although the multiplicity fluctuations are some- 
what overestimated for Mx larger than about 20 GeV. 

The present analysis adds new support for the conclusion, derived from studies of event 



shapes H 12] and from a study of energy flow and single particle momentum spectra || in 



large-rapidity-gap events in HI, that gluons play a prominent role in deep-inelastic diffraction. 
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Figure 1: Multiplicity moments (n), D and R2 in full phase space (a-c) and in single hemispheres 
(d-f) for charged hadrons as a function of Mx (HI), Mx (meson diffraction), W (EMC), yfs 
(e + e~), respectively. For clarity, HI data points in single hemispheres are slightly shifted in the 
horizontal direction with respect to their true positions. Also shown are predictions of several 
Monte Carlo models (see text). The Monte Carlo curves in forward and backward hemispheres 
are symmetric and their average is plotted. 
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Figure 2: The multiplicity distribution in KNO form for three intervals in Mj (DIS data and 
Monte Carlo) and at ^fs = (Mx) (JETSET e + e~), in full phase space for positive and negative 
particles separately (a-c), and for all charges in single hemispheres (d-f). The error bars show 
statistical errors only. Also shown are predictions of several Monte Carlo models (see text). The 
Monte Carlo curves are charge and forward-backward symmetric and their average is plotted. 
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Figure 3: Charged particle rapidity spectra for three intervals in Mx (HI), at W = (Mx) {fiN) 
and at y/s = (M x ) (JETSET e + e"). The (W) values for EMC and E665 differ slightly from the 
ones indicated for HI and are equal to 5.2, 11 and 19 GeV and 11.4 and 23.6 GeV, respectively. 
Also shown are predictions of several Monte Carlo models (see text). 
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Figure 4: (a) Central region charged particle density (—0.5 < y < 0.5); (b) The parameter p 
(HI) and b (others), which reflects the correlation between the number of particles in the forward 
and backward hemispheres, as a function of Mx (HI), W = (Mx) (pp), Mx (hadron diffraction), 
t/s (non-diffractive hadron-hadron) , yfs (JETSET e + e~). Also shown are predictions of several 
DIS Monte Carlo models (see text). 
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